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Genetic diversity and antimicrobial resistance 
in Streptococcus agalactiae strains recovered 
from female carriers in the Bucharest area 
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For the first time, we used multilocus sequence typing (MLST) to understand how Romanian group B strepto- 
coccus (GBS) strains fit into the global GBS population structure. Colonising isolates recovered from adult human 
females were tested for antibiotic resistance, were molecularly serotyped based on the capsular polysaccharide 
synthesis (bps) gene cluster and further characterised using a set of molecular markers (surface protein genes, 
pilus-encoded islands and mobile genetic elements inserted in the scpB-lmb intergenic region). Pulsed-field gel 
electrophoresis was used to complement the MLST clonal distribution pattern of selected strains. Among the 55 
strains assigned to six cps types (la, lb, II-V), 18 sequence types (STs) were identified by MLST. Five STs represented 
new entries to the MLST database. The prevalent STs were ST-1, ST-17, ST- 19 and ST-28. Twenty molecular marker 
profiles were identified. The most common profiles (rib+GBSil+PL-1, nb+GBSil+PL-1, PI-2b and alp2/3+P/-7, PI- 
2a) were associated with the cps LLI/ST-17 and cps V/ST-1 strains. A cluster of fluoroquinolone-resistant strains was 
detected among the cps V/ST-1 9 members; these strains shared alpl and IS154S and carried PI-1, PL-2a or both. Our 
results support the usefulness of implementing an integrated genotyping system at the reference laboratory level to 
obtain the reliable data required to make comparisons between countries. 
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Despite prophylactic perinatal antibiotic treatment, 
which has significantly diminished the burden of early 
onset disease, Streptococcus agalactiae [group B strep- 
tococcus (GBS)] remains a major infectious threat in 
neonates. Consequently, the development of a vaccine 
against GBS is a priority in health care research. Ef- 
forts to develop a vaccine capable of conferring broad 
coverage are also justified by the new challenges arising 
in GBS epidemiology, such as the diversification of the 
GBS serotypes (Slotved et al. 2007), a shift in serotype 
prevalence (Diedrick et al. 2010), capsular switching 
among strains (Bellais et al. 2012), a rising incidence of 
invasive disease in adults (Edwards & Baker 2005) and 
a growing rate of reduced GBS susceptibility to antimi- 
crobials (Castor et al. 2008). In this context, the monitor- 
ing of circulating strains and their characteristics using 
laboratory surveillance is pivotal for current and future 
vaccine development strategies. 

As knowledge about the GBS genome has accumu- 
lated, genotypic systems have been developed and im- 
proved to investigate the population structure of this 
organism and to distinguish differences between strains 
isolated from different sources. Pulsed-field gel electro- 
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phoresis (PFGE) and multilocus sequence typing (MLST) 
systems established for GBS have proven to be sufficient- 
ly discriminatory for epidemiological studies (Fasola et 
al. 1993, Jones et al. 2003). MLST was subsequently asso- 
ciated with a public website and a S. agalactiae database 
and has become a useful method for unifying global epi- 
demiology data on asymptomatic carriage of and infection 
by this pathogen. MLST revealed distinct GBS sequence 
types (STs) that were strongly associated with invasive 
disease cases. Of note, the ST-17 lineage has been associ- 
ated with neonatal infection, which facilitated the devel- 
opment of a rapid-screening molecular assay for routine 
diagnostic practice (Lamy et al. 2006). Population-based 
surveillance studies in different regions of the world have 
shown that the distribution of the GBS genetic lineages 
varies by population, which could limit the effectiveness 
of a universal GBS vaccine (Marchaim et al. 2006, Bohn- 
sack et al. 2008, Manning et al. 2009). 

Few molecular epidemiology studies of S. agalac- 
tiae have been undertaken in our region, mainly due 
to the high costs and the logistical difficulties that are 
encountered. Consequently, the paucity of information 
has hindered us from discerning the distinctive genomic 
traits of the local isolates and understanding how these 
isolates fit into the global population structure of this or- 
ganism. This study is the first to use MLST genotyping 
to investigate the clonal relationships and lineage distri- 
butions of GBS isolates recovered from female carriers 
residing in the Bucharest metropolitan area. In addition 
to MLST, antibiotic resistance and the expression of a set 
of molecular markers associated with the capsular poly- 
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saccharide synthesis (cps) gene cluster, virulence-asso- 
ciated surface protein antigen genes and mobile genetic 
elements (MGEs) were comprehensively characterised. 
The distributions of pathogenicity islands (PI)-l and 
PI-2 (variants PI-2a and PI-2b), which carry the genes 
required for the biosynthesis of pili (Telford et al. 2006), 
were also investigated because of the importance of pili 
in bacterial pathogenesis and vaccine development strat- 
egies (Lauer et al. 2005, Margarit et al. 2009). 

SUBJECTS, MATERIALS AND METHODS 

Bacterial strains - A total of 55 non-invasive GBS 
strains from the Bucharest metropolitan area were in- 
cluded in this study. These strains were recovered from 
vaginal swab specimens (1 isolate per specimen) of 55 ep- 
idemiologically unrelated, nonpregnant adult women and 
the samples were received by a large, private clinical labo- 
ratory in the community. Forty-three strains collected in 
2011 were fully characterised in this study. Another set of 
12 strains, collected in 2010, were only subjected to MLST 
analysis and pilus island-associated gene assessment, as 
their antimicrobial resistance, capsular genotype, surface 
protein encoding genes and MGEs were defined in a pre- 
vious study (Usein et al. 2012). Four of these strains that 
were identified as ST-17 by the rapid polymerase chain 
reaction (PCR) assay were validated using the MLST 
method to rule out any genotypic discrepancy. 

Template DNA - Template DNA was prepared from 
GBS cultures using a Chelex-based boiling procedure. GBS 
cells were cultured overnight on Columbia agar base (Ox- 
oid, Basingstoke, UK) supplemented with 5% defibrinated 
sheep blood. Then, the cultures were suspended in 200 uL 
of 10% Chelex 100 resin (BioRad) dissolved in TE buffer 
(10 mM Tris-HCl, 1 mM EDTA, pH 8.0); the cultures were 
boiled for 10 min and centrifuged for 5 min at 13,000 rpm. 
After centrifugation, 100 uL of the supernatant, which 
contained total bacterial DNA, was removed and stored at 
-20°C until further use in the molecular assays. 

Antimicrobial susceptibility - The isolates were test- 
ed for susceptibility to penicillin (10 U), erythromycin 
(15 ug), clindamycin (2 itg), tetracycline (30 Ltg), levo- 
floxacin (5 Ltg) and moxifloxacin (5 Ltg) using the disk 
diffusion method; the European Committee on Antimi- 
crobial Susceptibility Testing susceptibility breakpoints 
were used for all of the antibiotics (Oxoid) (eucast.org/ 
clinical breakpoints). Isolates sensitive to clindamycin, 
but resistant to erythromycin, were further tested using 
the double-disk diffusion or D-zone test to determine 
whether resistance to clindamycin could be induced 
and these isolates were classified as having either the 
inducible macrolide, lincosamide and streptogramin B 
(iMLS B ) antibiotic resistance phenotype or the M pheno- 
type. Isolates resistant to clindamycin and erythromycin 
were confirmed as the constitutive MLS B phenotype. 

PCR-based protocols using previously described 
primers were performed to detect the presence of the 
macrolide and tetracycline-resistance genes: ermA, 
ermB, mefA, tetM, tetO and tetL (Usein et al. 2012). 

The gyrA and parC gene fragments, including the 
quinolone resistance-determining regions (QRDRs), 



were PCR amplified from the GBS isolates that were 
resistant to levofloxacin and moxifloxacin using previ- 
ously published primers (Kawamura et al. 2003). The 
474 and 586 bp amplification products of gyrA and 
parC, respectively, were sequenced in both directions 
using the same primers. DNA sequencing reactions were 
performed with the Big Dye Terminator v.3.1 Cycle Se- 
quencing kit (Applied Biosystems, Foster City, Califor- 
nia, USA) and were analysed in an ABI 3130 sequencer 
(Applied Biosystems). The nucleotide and deduced pro- 
tein sequences were analysed using publicly available 
software (BioEdit) and compared with the S. agalactiae 
A909 gyrA and parC sequences, which are available in 
the GenBank database (accession NC 007432.1). 

Molecular typing - MLST - The internal fragments of 
seven housekeeping genes (adhP, pheS, atr, glnA, sdhA, 
glcK, tki) were amplified for MLST, according to the GBS 
scheme described by Jones et al. (2003). All PCR products 
were purified using a commercial kit (Wizard SV Gel and 
PCR Clean-Up System, Promega, Madison, Wisconsin, 
USA) according to the kit instructions and sequenced on 
both strands. Sequencing was performed with the ABI 
Prism BigDye Terminator v.3.1 Cycler Sequencing kit 
(Applied Biosystems). Before subsequent capillary elec- 
trophoresis on an ABI 3130 Genetic Analyzer (Applied 
Biosystems), the sequencing products were purified using 
the DyeEx 2.0 Spin kit (Qiagen GmbH, Hilden, Germa- 
ny). Raw sequence traces were reviewed by visual inspec- 
tion using the BioEdit software. The consensus sequences 
were submitted to the S. agalactiae database (pubmlst. 
org/sagalactiae/) to obtain an allelic profile and to assign 
a ST. The MLST data generated by the study isolates were 
further analysed using the eBURST software program, 
which can be found at pubmlst.org/analysis/. The default 
eBURST setting identified groups of related STs using 
the most stringent (conservative) definition, such that all 
members assigned to the same group shared identical al- 
leles at six of the seven loci with at least one other member 
of the group. The study STs were compared with all STs 
in the S. agalactiae MLST database. A phylogenetic tree 
was drawn from the allelic profile data using the Phylo- 
dendron software (pubmlst.org/analysis/), which employs 
an unweighted pair group method with an arithmetic av- 
erages (UPGMA) algorithm. 

Sequencing of the modified glcK allele - In addition 
to the glcK amplification and sequencing primers from 
pubmlst.org/sagalactiae/, specific primers were designed 
to sequence the insert detected within the glcK amplicon 
of one strain using the Primo primer design tool (mob- 
yle.pasteur.fr); the primers were based on the sequence 
of a previously published glcK allele available in Gen- 
Bank (accession EF990366) (Martins et al. 2007). These 
primers are listed in Table I. 

PFGE - The strains were compared by PFGE of Smal- 
digested genomic DNA, using a protocol described else- 
where (Usein et al. 2012), on a CHEF-Mapper apparatus 
(Bio-Rad Laboratories, Hercules, CA, USA). Cluster 
analysis of the PFGE fragment patterns was performed 
using the Dice coefficient at 1.5% tolerance and 1.5% 
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TABLE I 

Primers designed in this study 
for sequencing the insert of glcK allele 



Nucleotide sequence 
Primer (5'->3') 



glcK 


„362F 


CAGGTGGAGAAATTGGGC 


glcK_716F 


GATGGCACCTAAGGTATTCAC 


glcK 


.1079F 


TCAGGATTGTCAAGGCTACTC 


glcK 


1518R 


AGGAGAAGCTTTTCGTGATAG 


glcK 


_1882R 


TCAATTCTTTAGCCTCTAATGC 


glcK 


_2236R 


CCTTTCTTTGGATGACGACG 


glcK 


_2769R 


TGTCGCTGATGCAACTG 



optimisation and a dendrogram was constructed using 
the UPGMA algorithm and the Bionumerics software 
(Applied Maths, Sint-Martens-Latem, Begium). The re- 
sults were interpreted according to the criteria proposed 
by van Belkum et al. (2007). Strains with PFGE profiles 
that differed by four or fewer bands were assigned to 
subtypes of the same PFGE type (pulsotype) and isolates 
with patterns that differed by five or more bands were 
assigned to distinct types. 

cps genotype - A previously published, 19-plex PCR 
assay for the identification of the la, lb and II-LX GBS 
capsular types was used, with slight modifications (Im- 
peri et al. 2010). PCRs were performed in a final vol- 
ume of 25 (O.L. Each PCR reaction mixture contained 5 
uL of bacterial lysate, lx Qiagen Multiplex PCR Master 
Mix (Qiagen, Hilden, Germany) and the 19-primer pool 
(primers 2-15, 17-19 at 6 pmoles each and primers 1 and 
16 at 10 pmoles each). Amplification was carried out in a 
GeneAmp PCR System 2700 (Applied Biosystems) and 
the following thermocycling conditions were used: an 
initial denaturation of 95°C for 15 min (to activate Hot- 
StarTaq DNA polymerase) followed by 15 cycles of 94°C 
for 30 s, 54°C for 90 s and 72°C for 90 s, then 25 cycles of 
94°C for 30 s, 56°C for 90 s and 72°C for 90 s and a final 
extension at 72°C for 10 min. 

Molecular traits - Surface protein coding gene pro- 
file and MGE detection - Previously described primers 
and protocols were used to identify the presence of the 
bac, bca, rib, alpl, alp2/3 and alp4 genes that encode 
C beta, C alpha, Rib and the alpha-like proteins Alpl, 
Alp2, Alp3, Alp4, respectively, as well as the insertion 
of the MGEs IS1548 and GBSil within the scpB-lmb ge- 
netic region (Usein et al. 2012). 

Pilus-coding gene profile - Previously published 
primers were used to screen for the presence of PI-1, PI- 
2a and PI-2b (Margarit et al. 2009). The PI-1 backbone 
protein-coding gene (1,798 bp expected amplicon), the 
PI-2a ancillary protein-coding gene (-2,840 bp expected 
amplicon) and the PI-2b backbone-coding gene (1,700 bp 
expected amplicon) were selected as targets for the iden- 
tification of the PI-1, PI-2a and PI-2b-like regions. 



The PCRs had a final reaction volume of 50 uE, 
which contained 5 uE of bacterial lysate, lx High Fi- 
delity PCR Master (Roche, Mannheim, Germany) and 
15 pmoles each of the forward and reverse primers. The 
PCRs were run on a GeneAmp PCR System 2700 (Ap- 
plied Biosystems), using the following thermal cycling 
conditions: an initial denaturation of 94°C for 3 min fol- 
lowed by 9 cycles of 94°C for 20 s, 55°C for 60 s and 
72°C for 90 s, then 24 cycles of 94°C for 20 s, 55°C for 
30 s and 72°C for 90 s plus 5 s for each successive cycle 
and a final extension at 72°C for 7 min. 

RESULTS 

MLST analysis - The 55 strains examined in this 
study were resolved into 18 unique STs. The number of 
allelic profiles for each of the seven housekeeping genes 
of the GBS MLST scheme ranged from three (pheS) to 
11 alleles (atr). We identified four atr alleles (atr-63, atr- 
64, atr-65 and atr-66) and five STs (ST-587 and ST-592 
through 595), which were newly assigned by the data- 
base curator. Four of them were variants related to ST-19 
(ST-594, ST-595), ST-22 (ST-592) and ST-297 (ST-593) 
with new atr alleles. ST-587 had a new combination of 
previously known alleles. The sequences of all the new 
alleles and the composition of the newly identified STs 
are available from pubmlst.org/sagalactiae/. 

MLST genotyping revealed that six STs were repre- 
sented by at least three strains each: ST-1 (8 strains), ST- 8 
(3 strains), ST-17 (13 strains), ST-19 (9 strains), ST-23 (3 
strains) and ST-28 (6 strains). ST-297 was detected twice 
and 11 STs (ST-10, ST-12, ST-26, ST-49, ST-103, ST-335, 
ST-578 and ST-592 through ST-595) occurred only once. 
The presumed ST-17 status of the four strains from the 
previous study was confirmed by MLST. 

The Romanian strains' MLST profiles were includ- 
ed in four of the 12 groups defined by eBURST when 
compared with all known allelic profiles from the GBS 
MLST database, downloaded on 11 July 2013 (Fig. 1). 
Fifteen STs (ST-1, ST-8, ST-10, ST-12, ST-17, ST-19, ST- 
28, ST-49, ST-297, ST-335, ST-587 and ST-592 through 
ST-595) were categorised into group 1, the largest and 
most complex eBURST group. These STs clustered into 
multiple clonal complexes (CCs): CC19 (ST-19, ST-28, 
ST-335, ST-587, ST-594 and ST-595), CC1 (ST-1, ST-49, 
ST-297 and ST-593), CC10 (ST-10, ST-8 and ST-12), CC17 
(ST-17) and CC67 (ST-592). The remaining three STs 
(ST-23, ST-103 and ST-26) were distributed into groups 
2, 3 and 6, respectively, representing CC23, CC26 and 
CC103. When the autochthonous strains were analysed 
as an independent group, setting the group definition of 
seven shared alleles to zero, ST-17, ST-23, ST-26, ST-103 
and ST-592 appeared as individual, unlinked STs (sin- 
gletons), along with three clusters corresponding to the 
above CC19, CC1 and CC10 groups. 

Sequence of the modified glcK allele - For one of 
the strains, the internal fragment of the glcK allele gave 
a larger-than-expected PCR product, estimated at ap- 
proximately 3 kb instead of the predicted 0.459 kb. DNA 
sequencing of the amplicon revealed that the increase in 
band size was due to a 2,313 bp insertion between nu- 
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eBURST 
group 



. ST-l 
. ST^9 
. ST-297 

- ST-593 
. ST-19 

. ST-28 
. ST-587 
. ST-335 
. ST-594 
. ST-595 

- ST- 103 
. ST-8 

. ST-10 
. ST-l 2 
. ST-26 
.ST-17 
. ST-592 
. ST-23 



No. of isolates 



Pilus gene profile 
(no. of isolates) 
PI-l+PI-2a (8) 
Pl-l+PI-2a(l) 
PI-1 (1); PI-l+PI-2b(l) 
PI-l+PI-2b(l) 

PI-1 (1); PI-2a (1); PI-l+PI-2a (7) 

PI-2a(3);PI-l+PI-2a (3) 

PI-l+PI-2a(l) 

PI-l+PI-2a(l) 

PI-l+PI-2a(l) 

PI-l+PI-2a(l) 

PI-2b(l) 

PI-l+PI-2a (3) 

PI-l+PI-2a(l) 

PI-l+PI-2a(l) 

PI-2a(l) 

PI-1 (1); PI-2b (6); PI-l+PI-2b (6) 
PI-2a (1) 

PI-2a (2); PI-l+PI-2a (1) 



Fig. 1 : unweighted pair group method with an arithmetic averages tree generated from the multilocus sequence typing allelic profiles of the 55 isolates 
in relation to the eBURST groups and pilus islands distribution. 



cleotides 405 and 406 of the gkcK fragment. A BLAST 
search identified an open reading frame (ORF) within this 
sequence, with 98% identity to a group II intron reverse 
transcriptase. The newly determined nucleotide sequence 
was deposited in GenBank (accession HF679144). The 
incorporated sequence was removed for MLST analysis, 
leaving an intact glcK fragment identified as allele 2. The 
complete allelic profile assigned the isolate to ST-l. 

Presence of pilus islands and correlations with STs - 
Genes associated with at least one of the three GBS pilus 
islands were detected in each of the 55 strains investigat- 
ed. More precisely, PI-1 was detected in 40 strains, PI-2a 
in 37 strains and PI-2b in 15 strains. Nineteen strains 
displayed only one type of pilus island and 36 strains 
carried a two-pilus island combination. Consequently, 
the PCR results were considered predictive for the pres- 
ence of the following profiles: PI-1 (3 strains), PI-2a (8 
strains), PI-2b (7 strains), PI-l+PI-2a (29 strains) and PI- 
l+PI-2b (8 strains). The concurrent presence of PI -2a and 
PI-2b was not detected. For the STs represented by more 
than three strains, certain correlations were observed 
between the presence of a particular pilus island profile 
and the ST status. ST-l strains exclusively displayed PI- 
l+PI-2a profiles. With the exception of an ST-19 strain 
that carried PI-1 alone, the ST-19 and ST-28 strains con- 
tained PI-2a alone or in association with PI-1. All ST- 
17 strains but one possessed PI-2b and this island was 
detected in only three other STs, namely ST-103, ST-297 
and its single locus variant (SLV), ST-593. 

In terms of the pilus island distribution across the sero- 
types, the PI-l+PI-2a profile was present in strains display- 
ing cps genotypes V (17 strains), III (5 strains), II (5 strains) 
and lb (2 strains). PI-2b, with or without PI-1, was identi- 
fied in cps III (12 strains), IV (2 strains) and la (1 strain). 

Identification of cps genotypes, surface protein en- 
coding genes andMGEs - Molecular serotyping success- 



fully assigned all 43 strains uncovered in this study to 
six cps types: la (3 strains), lb (1 strain), II (7 strains), III 
(13 strains), IV (2 strains) and V (17 strains). 

PCR amplification of the surface protein encoding 
genes revealed that 41 strains were positive for one of 
the following genes: bca (4 strains), alpl (12 strains), 
alp2/3 (7 strains) and rib (18 strains). Three of the bca- 
positive strains concurrently carried the CP protein gene, 
bac. Two strains were PCR-negative for all of the sur- 
face protein-encoding genes investigated. Twenty-nine 
strains (67%) harboured either the group II intron GBSil 
(20 strains) or the insertion sequence IS1548 (9 strains) 
within the scpB-lmb genetic region. 

The 12 strains that were obtained from a previous 
study were assigned to the following cps types: type lb 
(1 strain), type II (2 strains), type III (5 strains), type 
IV (1 strain) and type V (3 strains). These strains har- 
boured the following surface protein encoding genes: 
bca (1 strain), alpl (2 strains) alp2/3 (2 strains) rib (7 
strains) and bac (1 strain). The bca and bac genes were 
associated. Four of the 12 strains contained neither of the 
two examined MGEs, whereas seven were positive for 
the intron GBSil and one carried IS1548. 

Table II presents the relationships between the se- 
quence types, cps genotypes and the combinations of 
various genetic markers (surface protein genes, MGEs, 
pilus islands) for all 55 strains. 

Antimicrobial resistance profiles, resistance geno- 
types and associations with cps type - All 43 isolates 
were susceptible to penicillin, but resistant to tetracy- 
cline either due to the presence of tetM alone (35 strains) 
or due to tetM in association with tetO (8 strains). Eryth- 
romycin resistance was observed in 19 strains. Nine of 
these strains displayed a constitutive clindamycin resis- 
tance phenotype due to ermB, while 10 strains expressed 
an inducible phenotype due to ermA alone (5 strains), 
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TABLE II 

Relationships between sequence types, capsular polysac- 



charide synthesis (cps) genotypes and the molecular marker 


profiles (i 


.e. surface protein genes, mobile genetic elements, 




pilus islands) for the 55 studied strains 






IVlUlCCUlal IllallvCI piUlllC 


Sequence 


(isolates) 


(isolates) 


type 


(n) 


(n) 


1 


V(8) 


alp2/3 + PI-1, PI-2a (7) 






rib + PI-1, PI-2a (1) 


8 


lb (1) 


bac, bca + PI-1, PI-2a (1) 




V(2) 


bac, bca + GBSil + PI-l,PI-2a (2) 


10 


II (1) 


alpl +GBSil + PI-1, PI-2 a (l) 


12 


lb (1) 


bac, bca + PI-1, PI-2a (1) 


17 


111(13) 


rib + GBSil + PI-1, PI-2b (6) 






rib + GBSil + PI-2b (6) 






rib + GBSil + PI-1 (1) 


19 


III (2) 


rib + 1S1548 + PI-1, PI-2a (2) 




V(7) 


alpl + 1S1548 + PI-1, PI-2a (3) 






alpl + GBSil + PI-1, PI-2a (2) 






alpl + 1S1548 + PI-1 (1) 






alpl + IS1548 + PI-2a (1) 


23 


la (2) 


alpl + PI-2a (2) 




HI (1) 


alp2/3 + PI-1, PI-2a (1) 


26 


V(l) 


GBSil + PI-2a (1) 


28 


II (6) 


rib + GBSil + PI-1, PI-2a (3) 






mfcXTDCil _l_ T>T "in f"2\ 

no + Oooii > ri-za {3) 


49 


V(l) 


alp2/3 + PI-1, PI-2a (1) 


103 


Ia(l) 


PI-2b (1) 


297 


IV (2) 


alpl + PI-l,PI-2b(l) 






alpl + PI-1 (1) 


335 


III (1) 


rib + 1S1548 + PI-1, PI-2a (1) 


578 


11(1) 


rib + GBSil + PI-1, PI-2a (1) 


592 


II (1) 


bca + 1S1548 + PI-2a (1) 


593 


IV (1) 


alpl + PI-1, PI-2b (1) 


594 


HI (1) 


rib + 1S1548 + PI-1, PI-2a (1) 


595 


V(l) 


alpl + GBSil + PI-1, PI-2a (1) 



ermA associated with mefA (4 strains) or ermB (1 strain). 
Five of the strains with inducible resistance to clin- 
damycin were concurrently resistant to levofloxacin and 
moxifloxacin (Fig. 2). These strains displayed the cps V 
genotype and the parC and gyrA sequences, including 
the QRDRs responsible for fluoroquinolone (FQ) resis- 
tance, showed an identical pattern of nucleotide muta- 
tions. Specifically, the strains presented single substitu- 
tions in GyrA (Ser81Leu) and ParC (Ser79Phe), as well 
as identical silent mutations that were found exclusively 
in parC (codons 142 and 145). 

PFGE analysis ofFQ-resistant (FQR) isolates - FQR 
strains were also subjected to PFGE analysis. All strains 
yielded interpretable PFGE profiles following Smal 
digestion and these profiles were compared using an 
UPGMA-based dendrogram (Fig. 2). The FQR cluster 



extended to 72% similarity, corresponding to a total of 
four DNA fragments or bands that differed. Consequent- 
ly, these strains were assigned to the same pulsotype ac- 
cording to the recommendations of van Belkum et al. 
(2007). Four subtypes were identified within this pulso- 
type, with 86% similarity and two bands that differed. 
Two strains were considered indistinguishable and were 
assigned to the same subtype. 

DISCUSSION 

This is the first Romanian study to investigate the 
diversity of the GBS strains colonising the human fe- 
male population using the MLST approach linked to the 
analysis of significant genetic markers for GBS descrip- 
tion. This comprehensive genotyping system generated 
a highly informative characterisation of autochthonous 
strains and allowed for objective comparisons with 
strains from other geographic regions. Overall, the GBS 
population sampled in this study comprised 18 ST pro- 
files, six capsular genotypes and 20 molecular profiles 
defined by the combination of three sets of markers (6 
surface protein genes, 2 MGEs and 3 pilus islands). 

Most STs identified in this study have been previously 
described; however, despite the small number of strains, 
we also identified five STs that were not present in the 
MLST database. The most prevalent STs in our data set 
were ST-1, ST-17, ST-19 and ST-28. These results were 
partially in agreement with a study of a diverse, globally 
derived GBS collection of 67 strains isolated from adult 
carriers (Jones et al. 2003). In that study, the ST-1, ST-17 
and ST-19 strains also occurred frequently, but the ST-28 
clone was less prevalent. By contrast, a higher propor- 
tion of ST-23 members was reported, a finding that was 
also reported in a study of colonising strains isolated 
from young nonpregnant American women (Manning et 
al. 2008). In the latter study, ST-22 was the second-most 
prevalent clone identified after ST-1, followed by ST-23 
and ST-19. ST-22, a clone that also frequently occurred 
in pregnant women from Israel (Bisharat et al. 2005), 
was not identified among the Romanian samples, except 
for one strain that was assigned to ST-592, an ST-22 SLV. 
ST-19 and ST-28 belong to the same CC19 group. ST- 
28 strains were more prevalent within CC19 in a study 
focusing on strains isolated from maternal carriage in 
two different African cities, Dakar and Bangui (Brochet 
et al. 2009). In the Romanian collection, as well as in 
others from Europe and the United States, ST-19 was de- 
scribed as the most prevalent ST included in this CC. 

The ST-17 clone, which is significantly associated 
with invasive infections, was strongly represented in the 
present study. A frequent occurrence of lineage ST-17 
colonisation has also been reported for pregnant East 
African women, German adults and Jews and non-Jews 
who emigrated from the former Soviet republics (Mar- 
chaim et al. 2006, Eickel et al. 2009, Huber et al. 2011). 

For STs with many isolates, such as ST-1, ST-17, ST- 
23 and ST-28, we found several associations between the 
MLST profiles and the molecular markers screened in 
this study. Among this collection's strains, the cps II, 
III and V types predominated. In agreement with other 
studies, although a capsular serotype was generally not 
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restricted to a specific ST, we observed the tendency for 
a single serotype to prevail within individual STs (Bro- 
chet et al. 2006, Bohnsack et al. 2008). In this study, the 
exclusive associations ST-l/cps V, ST-lllcps III and ST- 
28/cps II should be mentioned. In contrast with other 
studies that reported a strong association of ST-19 with 
serotype III, in our collection, the members of this clone 
were more frequently associated with serotype V (Mar- 
chaim et al. 2006, Manning et al. 2009). 

Also of note, one of the ST-l/cps V strains in this 
study carried a putative group II intron inserted within 
the glcK internal fragment, which was discovered in the 
GBS MLST scheme. While searching the literature for 
similar findings, we found three other studies that men- 
tioned isolates harbouring the same type of insertion, 
each located exclusively in a glcK allele in the same posi- 
tion (Bisharat et al. 2004, Brochet et al. 2006, Martins et 
al. 2007); however, the sequence data that were available 
for comparisons were limited. The previously described 
isolates were mainly of bovine origin and belonged to 
other lineages (e.g., ST-61). The only reported strain of 
human origin belonged to serotype la and the ST-248 
lineage (Brochet et al. 2006). 

Correlations were observed between the cps types 
and the surface protein genes. The cps III and II type 
strains were closely associated with rib and the cps lb 
strain was associated with bac and bca, which was con- 
sistent with other reports (Kong et al. 2002, Martins et 
al. 2011). The association of serotype V with alp3 that 
was evident in other studies (Kong et al. 2002, Martins 
et al. 2011) was not encountered in our collection, where 
a similar number of strains displaying the cps V geno- 
type harboured either alp2/3 or alpl. 

Previous studies revealed some associations between 
certain STs and surface protein genes. In addition to the 
expected ST-17/r/& association, we observed that the ST- 
2%lcps II Romanian strains also exclusively harboured 
rib. In contrast, the ST-28 strains of serotype II from a 
German colonising collection carried either bca or alp2/3 
(Eickel et al. 2009). Two serotypes were detected in the 
ST-8 strains, but both serotypes carried the same asso- 
ciation with bac and bca. This was not the case for ST-19 
and ST-23, both of which included serotypes that har- 



boured different surface protein encoding genes. Inter- 
estingly, the only representative of ST-26 from our col- 
lection lacked surface protein genes, which was similar 
to results found in a Polish study (Sadowy et al. 2010). 

Al Safadi et al. (2010) reported that most CC1 strains 
did not harbour either GBSil or IS1548 within the scpB- 
Imb intergenic region and mainly found GBSil in CC17 
strains of serotype III and CC19 strains of serotype II. 
Our results showed both similarities to and differenc- 
es from that study. We also noticed that the CC1 iso- 
lates (ST-1, ST-49, ST-297 and ST-593) lacked MGEs, 
whereas all of the ST-17 isolates exclusively contained 
the GBSil intron. Within the CC19 cluster, the ST-28/c/tt 
II strains displayed GBSil, whereas most ST-19 isolates 
(78%) harboured IS1548. The preponderance of the ST- 
19/IS1548 isolates with the cps V genotype was mainly 
due to a group of isolates with genetic homogeneity, 
according to PFGE typing, that acquired resistance to 
FQs. These strains were resistant to tetracycline and 
macrolides, with an iMLS B phenotype correlated with 
the presence of ermA and mefA. Double and even triple 
combinations of erm and/or mefA have been reported by 
other authors (Culebras et al. 2002, Fluegge et al. 2004), 
but these combinations were not found in our Romanian 
GBS isolates. FQ resistance was supported by key muta- 
tions detected in the QRDRs of gyrA and parC, which 
encode the GyrA and ParC subunits of the DNA gyrase 
and topoisomerase IV, respectively. The inferred altera- 
tions of these enzymes seemed to be sufficient to cause 
clinical resistance to levofloxacin and moxifloxacin 
(Tazi et al. 2008). The assumption of a spread of multi- 
drug resistant GBS isolates based on our study would be 
misguided, due to the small number of isolates studied; 
however, the presence of a predominant pulsotype might 
be indicative of the local circulation of genetically re- 
lated FQR isolates. Based only on the results of PFGE 
analysis, a Japanese study concluded that FQR isolates, 
mainly of serotype lb, might belong to a single clone 
that acquired FQ resistance and rapidly disseminated 
throughout Japan (Murayama et al. 2009). Further, a re- 
cent report from China revealed that a large percentage 
of FQR GBS was due to serotype III isolates of the ST-19 
lineage (Wang et al. 2013). 
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In this study, for the first time, we screened Roma- 
nian GBS strains for the presence of the genetic deter- 
minants of pilus-like structures, which are considered 
promising vaccine candidates (Margarit et al. 2009). All 
of the studied strains carried at least one pilus island and 
the most common pilus island pattern consisted of the 
concomitant presence of the PI-1 and PI-2a loci, as previ- 
ously reported for colonising strains from Europe or the 
United States of America (Margarit et al. 2009). Across 
our collection, this combination was identified in more 
than half of the studied strains (53%) belonging to 12 STs. 
For prevalent STs, such as ST-1 and ST- 19, the PI-l+PI- 
2a genotype prevailed. In a recent study, Martins et al. 
(2013) noted the dominance of the PI-l+PI-2a combination 
in several worldwide-disseminated lineages, as well as its 
nearly ubiquitous presence in isolates belonging to CC19 
and its association with maternal colonisation. While 69% 
of the local strains carried PI-2a alone or in combination 
with PI-1, the PI-2b and PI-l+PI-2b profiles were detected 
in fewer isolates (27%), mainly in those of the serotype III 
and ST- 17 lineage. The significant association of PI-2b, 
alone or in combination with PI-1, with this particular lin- 
eage has been previously reported (Madzivhandila et al. 
2013) and it was suggested that the acquisition of this vari- 
ant might be partly responsible for this lineage's particular 
virulence and tropism (Martins et al. 2013). 

The data acquired in this study shed new light on 
the need for a more rigorous characterisation of the GBS 
isolates circulating in our area and on the value of ex- 
panding the laboratory typing toolset by including meth- 
ods such as MLST, which are helpful for both local and 
global GBS epidemiology. 
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